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ABSTRACT 

Lens models appropriate for representing cusped galaxies and clusters are developed. 
The analogue of the odd number theorem for cusped density distributions is given. 
Density cusps are classified into strong, isothermal or weak, according to their lens- 
ing properties. Strong cusps cause multiple imaging for any source position, whereas 
isothermal and weak cusps give rise to only one image for distant sources. Isothermal 
cusps always possess a pseudo-caustic. When the source crosses the pseudo-caustic, 
the number of images changes by unity. 

Two families of cusped galaxy and cluster models are examined in detail. The 
double power-law family has an inner cusp, followed by a transition region and an 
outer envelope. One member of this family - the isothermal double power-law model - 
possesses an exceedingly scarce property, namely, the lens equation is exactly solvable 
for any source position. This means that the magnifications, the time delay and the 
lensing cross-sections are all readily available. The model has a three dimensional 
density that is cusped like r~ 2 at small radii and falls off like r -4 asymptotically. 
Thus, it provides a very useful representation of the lensing properties of a galaxy or 
cluster of finite total mass with a flat rotation curve. The second set of models studied 
is the single power-law family. These are single density cusps of infinite extent. The 
properties of the critical curves and caustics and the behaviour of the lenses in the 
presence of external shear are all discussed in some detail. 

Key words: gravitational lensing - galaxies: structure 
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1 INTRODUCTION 

The aim of this paper is to examine the properties of lens 
models with density cusps. High-resolution imaging of the 
the nuclei of early-type galaxies by the Hubble Space Tele- 
scope (e.g., Lauer et al. 1995) have provided abundant ex- 
amples of cusps. In almost all early-type galaxies, the loga- 
rithmic gradient of the surface brightness with respect to 
projected radius is constant right down to the very cen- 
tre. Galaxy cores are hardly ever observed. This observa- 
tional evidence for the ubiquity of central cusps is reinforced 
by theoretical suggestions of dynamical processes by which 
cusps can grow in the basins of potential wells (e.g., Bahcall 
& Wolf 1976; Faber et al. 1996; Evans & Collett 1997). Fur- 
ther, the influential numerical simulations by Navarro, Frenk 
& White (1996) of galaxy formation in hierachical clustering 
cosmogonies found striking evidence for a universal cusped 
density law for dark matter haloes and clusters. 

This change in our understanding of the structure of 
galaxies and clusters motivates our study of geometrically 
thin, centrally cusped gravitational lenses. Hitherto, the 
emphasis in gravitational lensing studies has largely been 
placed on non-singular models, such as softened isothermal 
spheres and ellipsoids or Plummer models (e.g., Blandford & 
Kochanek 1987; Kochanek 1996; Keeton & Kochanek 1997; 
Schneider, Ehlers & Falco 1992, chap. 8). For such non- 
singular lenses, it is well-known that the total number of 



images is odd and that the number of even parity images ex- 
ceeds the number of odd parity images by one (e.g., MacKen- 
zie 1985; Burke 1981; Schneider, Ehlers & Falco 1992, chap. 
5). This paper begins with an outline of general results on 
geometrically thin, centrally cusped lens models appropri- 
ate for galaxies and clusters. In Section 3, a new and simple 
family of axisymmetric cusped lenses, the double power-law 
lenses, is presented, together with an application to radial 
arcs in clusters. Section 4 studies one member of this fam- 
ily, the isothermal double-power lens, in detail. It is wor- 
thy of special scrutiny, as the lens equation is solvable for 
any source position. This property is possessed by only two 
other known circularly symmetric lenses (the Schwarzschild 
lens and the isothermal sphere). Section 5 examines a fam- 
ily of density cusps of infinite extent, with equipotentials 
that are similar concentric ellipses. The isothermal model in 
this family has already been extensively investigated (e.g., 
Kassiola & Kovner 1993). The remaining members have not 
been examined in detail as isolated lenses, although the sta- 
tistical properties of ensembles of these models have been 
studied in two classic papers (Blandford & Kochanek 1987; 
Kochanek & Blandford 1987). 

In all of the following sections, we use the notation and 
conventions of the admirable text-book of Schneider, Ehlers 
& Falco (1992). In particular, we use (2:1,12) as Cartesian 
coordinates in the lens plane and (t/1,2/2) as Cartesian coor- 
dinates in the source plane. The Poisson equation relating 
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Figure 1. Multiple image diagrams for strong, isothermal and weak cusps. Typical x — y curves are shown in dashed lines for a number 
of source positions. Multiple imaging takes place for the source position y\ if the line x — y\ intersects the deflection angle curve more 
than once. In the strong cusp, no matter how distant the source position, the observer sees at least two images. In the isothermal cusp, 
multiple imaging only takes place if the source is sufficiently close to the cusp. For sources far enough away, only one image is observed. 
In the weak cusp, three images are seen for source close to the cusp, one image for distant sources. The change-over in the number of 
images occurs as the source crosses the radial caustic, as shown. 



the convergence k and the deflection potential ip is scaled to 
read V 2 ip = 2k. The deflection angle O is just Vip. 

2 GENERAL THEOREMS 

To orient ourselves, we start with the simpler case of ax- 
isymmetric lenses with density cusps. Theorems for general 
cusped lenses are deduced in Section 2.2. 

2.1 Axisymmetric Lenses 

Let us consider a geometrically thin, axisymmetric lens with 
convergence n(R), where R = \x\. For an axisymmetric lens, 
the ray trace equation is one-dimensional, as all light rays 
emanating from the source and received by the observer lie 
in the common plane incorporating the observer, lens and 
source. Let us choose coordinates in the lens plane so that 
x = (x,0) and in the source plane so that y — (y,0). The 
surface density is assumed to be piecewise continuous ev- 
erywhere except at the central cusp R — 0. We suppose 
that the deflection angle vanishes asymptotically, or, equiv- 
alently, the convergence falls off at large radii faster than 
1/R, i.e., 

lim Rk(R) = 0. (2.1) 

Let the convergence diverge at the origin like 

lim k(R) = 0(i?~ 7 ). (2.2) 

R — >0 

Then: 

(1) If 1 < 7 < 2, multiple imaging occurs for any position 
of the source y, no matter how large. The model is said to 



have a strong cusp. 

(2) If 7 = 1, the cusp is isothermal. Multiple imaging occurs 
for some source positions. For sufficiently large values of y, 
there is only one image. For sufficiently small values of y, 
there are at least two images. 

(3) If < 7 < 1, multiple imaging occurs for some source 
positions. For sufficiently large values of y, there is only one 
image. For sufficiently small values of y, there are at least 
three images. The model is said to have a weak cusp. 

(4) If 7 < 0, the model is uncusped in projection. A nec- 
essary and sufficient condition for multiple imaging is that 
the central convergence k(0) must exceed unity (Schneider, 
Ehlers & Falco 1992, p. 236). We remark that models which 
are uncusped in projection may still possess a mild cusp in 
the three-dimensional luminosity density. 

These statements are readily proved. As the cusp is ap- 
proached, the deflection angle O = Vi/> behaves like 

lim 6 = 0(R- 1+1 ). (2.3) 

R — *0 

Near the cusp, the deflection angle vanishes if 7 < 1, tends 
to a constant value if 7 = 1 and diverges if 7 > 1. Now, 
the deflection angle is an odd function of x. So, this implies 
that O is discontinuous at x = for isothermal cusps, while 
O possesses a singularity at x — for strong cusps. The 
gradient of the deflection angle behaves like 

This diverges for all cusped models. 

Figure 1 shows multiple image diagrams for strong, 
isothermal and weak cusps in typical centrally condensed 
galaxy or cluster models. The graphs show the deflection 
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angle O, together with x — y lines for a number of values 
of the source position. The location of the images is given 
by the solutions of the lens equation, and therefore by the 
intersection of the x — y lines with the deflection angle. So, 
for strong cusps, it is evident from Fig. 1 that there are al- 
ways two images, irrespective of the distance of the source. 
Isothermal cusps have the curious property that the number 
of images can change by one, as the source position is moved 
outward across a curve that we shall call the pseudo-caustic. 
This is illustrated in the second panel of Fig. 1. This curios- 
ity has previously been noticed in analyses of the singular 
isothermal sphere and ellipsoid (Kovner 1987a; Kormann, 
Schneider & Bartelmann 1994; Schneider, Ehlers & Falco 
1992, p. 243). Weak cusps always exhibit multiple imaging 
for sufficiently small y. As the source is moved away, there 
comes a critical point when it crosses the radial caustic and 
the number of images diminishes by two. This is shown in 
the third panel of Fig. 1. 

Let us note that the observations of the centres of 
early-type galaxies indicate that strong, weak and isother- 
mal cusps all can occur. For example, table 2 of Faber et 
al. (1997) gives the cusp slope 7 for a sample of 61 el- 
liptical galaxies. The observationally fitted values of 7 lie 
in the range 0.0 < 7 < 1.21. NGC 1199 is an example 
of a galaxy with a strong cusp, while NGC 4467 has an 
isothermal cusp. Weak cusps are the most predominant in 
Faber et al.'s (1997) sample. The cosmological simulations 
of Navarro, Frenk & White (1996) suggest that dark matter 
haloes possess a universal density law with a logarithmi- 
cally singular surface density. Subsequently, numerical in- 
vestigations (Fukushige & Makino 1997; Moore, Governato, 
Quinn, Stadel & Lake 1997) together with theoretical ar- 
guments (Evans & Collett 1997) indicated that the origi- 
nal work had underestimated the severity of the cusp. The 
weight of the evidence now suggests that dark haloes have 
three-dimensional density cusps like p ~ r -4 ^ 3 . In projec- 
tion, this becomes k ~ r ~ 2 / 3 and so lies within the weak 
cusp regime. 

2.2 General Lenses and Index Theorems 

There is a general theorem for non-singular lenses that states 
that the total number of images is odd and that the num- 
ber of even parity images exceeds the number of odd parity 
images by one (MacKenzie 1985; Burke 1981; Fukuyama & 
Okamura 1997). The image type is classified by studying the 
Jacobian matrix of the Fermat potential cf>. If both the eigen- 
values are positive or both are negative, then the images are 
of type I and III respectively. Such images have even parity 
and correspond to minima or maxima of the Fermat poten- 
tial. If one of the eigenvalues is positive and one negative, 
the image is of type II. Such images have odd parity and cor- 
respond to saddle-points of the Fermat potential. The odd 
number theorem is usually stated in the form 



m — mi + nin = 1, 



(2.5) 



where m is the number of images of type I and so on. It 
follows from this that the total number of images (m + mi + 
mil) is necessarily odd. 

To set the notation, let us quickly re-derive the odd 
number theorem for non-singular lenses. Introducing 2 = 
xi + ix2, z — xi — 1X2, then the lens equation becomes just 
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Figure 2. The only possible two image topography (the limacon) 
of the Fermat surface for lenses with strong density cusps. The 
saddle point is marked by S and the + and — signs indicate regions 
where the Fermat surface is higher (H) and lower (L) than the 
saddle. An image of type II occurs at the saddle, an image of type 
I at the minimum. The density singularity is marked. 



(e.g., Witt 1990) 



dz 



= <fe = o. 



(2.6) 



Let G\ be a contour enclosing some point P in the complex 
z plane. The index of the map co = 4t{z) at P is 



IndOfeP) = - 



din a;, 



(2.7) 



where is the image of C z . So, the index is the number 
of times the image of P in the complex w-plane is encir- 
cled, as P in the complex z-plane is encircled once in the 
anti-clockwise direction. If C z does not enclose any zeroes 
or singular points, the index vanishes. It is straightforward 
to show (e.g., Schneider, Ehlers & Falco 1992, p. 174) that 
the index of an extremum of the Fermat potential is +1, 
whereas the index of a saddle point is —1. Provided the de- 
flection angle tends to zero at large radii, then it must be 
true asymptotically that 



<h=h\z-z e \ + o(\z\), 



(2.8) 



where z s — y\ +iy2 is the complex position of the source. So, 
the point at infinity is a pole and it has an index of unity. By 
shrinking the contour so that it becomes enmeshed with the 
singular points and zeroes, it follows that the sum of all the 
indices must be equal to unity. If the lens is non-singular, 
there are no poles and only the contributions of the zeroes 
must be taken into account, i.e., 

1 = Ind(<fe, oo) = Ind((fe, Pi) = m - n„ + n m . (2.9) 

i 

This is the odd number theorem and holds good only for 
non-singular lenses with convergence falling off faster than 
l/\z\ at large radii. 

What is the analogue of this theorem for general lenses 
with density cusps? A careful examination of the density 
singularity at the origin of the z-plane is required to establish 
its index. Near the origin, the lens equation can always be 
cast into the form 

<fe = \V -ZbI + ^ts-, (2.10) 



where A is a constant (for circular lenses) or a real func- 
tion of the phase (for non-axisymmetric lenses) and the 
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Table 1. Some simple potentials corresponding to the convergence of the double power-law family as given 
in eqs. (3.3) and (3.8). (Below, £ and n are natural numbers. The summands are empty if £ = or n = 0.) 



L = 1 N = 2n V(-R) 



L = l N = 2n + l ip{R) 



L = U N = 1 



L = 2£ + l N = 1 



8Ca 2 



2n + 1 

2Ca 2 

71+ 1 

8Ca 2 
2£+ 1 

2Ca 2 
^ + 1 



n 

log(l + y/Ri/° + l) - E 2fc 1 _ 1(1 + fil /V 1 / 



, , , i / pl/a \ fc-1/2 



sub-dominant terms have been neglected. Substituting z = 
eexp(i#), this becomes 

<fe = ijec X p(i6»)-z s |+Ae 1 - 7 cxp(ie). (2.11) 

Letting e — > so that the contour is wrapped tightly around 
the origin, we deduce that the first term is dominant in the 
case of weak cusps (0 < 7 < 1) and the second term is 
dominant in the case of strong cusps (1 < 7 < 2). In the 
limit, the vector field <fe is constant for weak cusps, so the 
the index of the origin is and the odd number theorem 
(2.5) still holds good. For strong cusps, a different situation 
obtains. The index of the origin is +1, as the vector field 
is radial in the limit e — * 0. Lens models with strong cusps 
have an even number of images, as 



m - nu + nm = 0. 



(2.12) 



The number of even parity images is equal to the number of 
odd parity images. 

The instance of isothermal cusps (7 = 1) is more subtle. 
Now both terms in (2.10) matter and the index of the origin 
depends on the position of the source z„. As e — > 0, the lens 
equation (2.11) becomes 

<h=±\z B \+Aexp(i0), (2.13) 

where A = ife which depends only on the phase for mod- 
els with isothermal cusps. If \z s \ < 2\ip^\, then the index of 
the origin is +1, otherwise the index of the origin is 0. This 
equation therefore defines a curve in the source plane - the 
pseudo-caustic. If the source is outside the pseudo-caustic, 
there is an odd number of images according to (2.5). Cross- 
ing the pseudo-caustic produces an additional, single image 
of type II initially at the origin. Once the source is within 
the pseudo-caustic, there is an even number of images ac- 
cording to (2.12). Note that the pseudo-caustic differs from 
a true caustic because the magnification of a point source 
is finite (not infinite) and because the number of images 
changes by one (not two). The analogues in the lens plane 
are the pseudo-critical curves. 

As is well-known (Blandford & Narayan 1986; Schnei- 
der, Ehlers & Falco 1992, p. 178), the topography of the Fer- 
mat surface for non-singular lenses can be classified in terms 
of the local structure at saddle points. For three image ge- 
ometries, the shape of the critical isochrone corresponds ei- 
ther to a lemniscate or a limacon. For five image geometries, 



there are six possible topographies given by joining lemnis- 
cates and limacons. These theorems still hold for lenses with 
weak density cusps. But, if the lens possesses a strong den- 
sity cusp, then the topography of the Fermat surface is more 
restricted. For two image geometries, the only possibility is 
the limacon illustrated in Fig. 2. One of the images occurs 
at the saddle point marked by S, one in the basin of the 
Fermat surface. The density cusp corresponds to the high- 
est point. For four image geometries, there are five possible 
topographies, given by joining at least one limacon to either 
a lemniscate or a limacon. 



3 THE DOUBLE POWER-LAW LENSES 

We now present a new family of axisymmetric lenses with 
central density cusps, the double power-law lenses. Section 
3.1 explains our strategy for locating the models with simple 
potential-convergence pairs. The following sections examine 
properties and applications of the models in detail. 



3.1 The Potential-Convergence Pair 

Let us consider the family of cusped lenses with surface den- 
sity S as a function of projected radius £ given by 

CI* 



£(£) = 



(3.1) 



Here, (a, j3, 7) are positive constants and £ c is a scale-length. 
By choosing the scale in the lens plane as £ c , then the conver- 
gence k is defined as (e.g., Schneider, Ehlers & Falco 1992, 
chap. 5) 



k{R) = 



S(i^ c 



^crit 



c 2 D s 



(3.2) 



Scrit 47rGDd-D ds 

where E cr i t is the critical surface density and D B ,Dd and 
Dds are the observer-source, the observer-lens and the lens- 
source distances respectively. For our model (3.1), this gives 

C 



k(R) 



Ky(Rl/<x + l)a(/3- 7 ) ' 



(3.3) 



At small radii, the convergence is cusped like i?~ 7 , whereas 
in the outer parts, it falls like R^ 13 . The parameter a con- 
trols the extent of the transition region between the cusp 
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Figure 3. The variation of the radii of the tangential and radial 
critical circles with total projected mass for various values of 7. 
Note that radial critical curves and caustics do not exist if 7 > 1, 
as is evident on referring back to the multiple imaging diagrams 
of Fig. 1. 



and the envelope. The lens models (3.3) are clearly inspired 
by the 'Nuker' profile introduced by Gebhardt et al. (1997) 
and Faber et al. (1997). These authors show that it provides 
a compact description of the surface brightness of early-type 
galaxies. The lens models (3.3) are also related to the three- 
dimensional elliptical galaxy models introduced by Hern- 
quist (1990) and studied by Zhao (1996). The normalisation 
constant C is related to the total mass M by 



M = 2Cn£ S crit aB(a[2 - 7] , a[/3 - 2] ) , 



(3.4) 



where B(x,y) — F(x)T(y)/T(x + y) is the complete Beta 
function. For convenience, we also define a dimensionless 
mass M such that 



M = 



M 



crit 



(3.5) 



For the total mass to converge, we must insist that [3 > 2 
and that 7 < 2. 

The solution of the two-dimensional Poisson equation 
is familiar from the theory of the Newtonian potential of 
circular cylinders as (e.g., Routh 1892; Ramsey 1940) 

f-E /> oc 

ip(R) = 2\og R I drrn{r) + 2 I dr r K (r) log r. (3.6) 

JO J R 

Of course, the first term is the potential of the circular cylin- 
ders interior to R, which attract as if their mass were con- 
centrated on axis. The second term is the potential of the 
circular cylinders exterior to R. Integrating (3.6) by parts, 
we obtain: 



i/j(R) =2\ogR / drrn(r) 



f 

Jo 



+ 2 



f> OC f> c 

JR r Jr 



(3.7) 



dr' r'«(r'). 



Substituting the convergence (3.3), we see that in general 
the potential is not elementary. Of course, there is an enor- 
mous advantage to working with models with at least sim- 
ple deflection angles and preferably simple potentials. So, 
let us first outline conditions under which these quantities 
are reducible to a finite number of elementary functions. It 
is helpful to introduce two constants L and N such that 



L = 2a[2 - 7] - 1, N = 2a[/3 - 2] - 1. 



(3.8) 



Making the substitution R = tan 2c * 8, then the deflection 
angle is 



e = 



dip _ 4Cq 

Jr ' 

2Ca 



R 



dd' sin L 8' cos N 8' 



(3.9) 



R 



B(L/2 + 1/2, N/2 + 1/2; sin 2 8). 



Here, B(x,y;z) is the incomplete Beta function (see e.g., 
Abramowitz & Stegun 1965; Gradshteyn & Rhyzik 1980). 
Amongst other instances, it reduces to a finite numer of 
elementary functions if either (1) L is an odd natural number 
and N is arbitrary, or vice versa, or (2) L is an even natural 
number (or zero) and N is any natural number (or zero), or 
vice versa. 

The potential can similarly be reduced to 



V> =4Ca log R I d8' sin L 8' cos N 8' 
'0 



(3.10) 



+8Ca 



AS' 



sin a' cos ( 



d(j> sin L (f> cos N < 



A consultation of the discussion in sections 2.51 and 2.52 
of Gradshteyn & Rhyzik (1980) reveals that the indefinite 
integration in (3.10) is elementary if L and N are natural 
numbers (or zero), at least one of which is odd. Let us note 
explicitly that if L and iV are both even natural numbers 
(or zero), the potential is not elementary as it ultimately de- 
pends on the transcendent integral (2.644.5) of Gradshteyn 
& Rhyzik (1980). 

Some simple examples of lens models with convergence 
(3.3) possessing elementary potentials are given in Table 1. 
Henceforth, we specialise to the case L = 1, N = 1 for defi- 
niteness. 

3.2 The L = 1, N = 1 Family 

Here, the potential-convergence pair is just 
M(2-7) 2M 



■ log(l + R 



2- 7 n 



(3.11) 



R-<(1 + R 2 ->) 2 ' T 2-7' 

If 1 < 7 < 2, the lens has a strong cusp. This becomes 
isothermal when 7=1. For < 7 < 1, the lens possesses 
a weak cusp. The instance 7 = is the uncusped Plum- 
mer model, discussed in detail in Schneider, Ehlers & Falco 
(1992, chap. 8). Models with 7 < are unrealistic, as the 
density is not centrally concentrated. 

The scaled deflection angle O at a displacement x in the 
lens plane is 

2Msign(x)|x| 1 " 7 



e = Vtp - 



i + |xi 2 -7 



(3.12) 



where |a;| — R. In the lens plane, the radius of the tangential 
critical circle Rt is given by 

Rt(l + Rt' 1 ) = 2MR 2 - 1 . (3.13) 

Generally, this implicit equation cannot be solved analyti- 
cally. It is exactly solvable when 7 = 1 to give 



Rt = ±Vl + 8M- i, 
when 7 = to give 
Rt = V2M-1. 



(3.14) 



(3.15) 
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The tangential critical circles are just the radii of the bright 
Einstein rings formed when the source is exactly coincident 
with the lens. The radial critical curves have radius R T given 
by 

(1 + Rff = 2M((1 - ■y)Rp - R 2 ' 2 ~<). (3.16) 



When 7 = 0, this is solvable to give 



R 2 



M (4 + M) — M — 1. 



(3.17) 



When 7 > 1, this equation has no positive root and there 
is no radial critical curve. As is evident from the multiple 
imaging diagrams of Fig. 1, isothermal and strong cusps do 
not have radial critical curves. Figure 3 shows the variation 
of the critical curves with total projected mass. The curves 
are labelled with the value of 7. The graph shows that as 
M — > 00, the tangential critical curves always tend to y/2M, 
which is the Einstein ring radius of a point mass (let us recall 
M is mass measured in units of the 27r^S cr it). 

The tangential caustic coincides with the origin in the 
source plane y t = 0. The radial caustic exists only for weak 
cusps and is 



Ur 



i? r (7 + 2R 2 -~<) 

7 - 1 + R 2 r< ' 



(3.18) 



The radial caustics separates the source domains (\y\ < y T ) 
generating three images from those (\y\ > y r ) generating 
just one. The magnification is 



2MR- 



1 



1 + R 2 —i 
2M((1 - 



j)R~ 



R 



2-2 7 \ 



(l + R 2 -<) 2 



(3.19) 



When I y \ < y T , the image at x > Rt is of type I, the image at 
— Rt < x < —R r is of type II and the image at — _R r < x < 
is of type III. The two components of the shear are 



7i 



72 



M{x 2 2 -x 2 1 ){ 1 + 2R 2 -~<) 
^Mx^^i + 2R 2 -"') 



(3.20). 



R 2 +t(l + R 2 - 



7^2 



3.3 An Application: Radial Arcs 

Radial arcs are images of galaxies distorted by foreground 
clusters. They are elongated in the radial direction, in dis- 
tinction to the more commonly occurring tangential arcs. 
Bartelmann (1996) described how, in a situation where both 
radial and tangential arcs are observed, the positions of the 
arcs can be used to determine the parameters of the model 
potential representing the lens. Bartelmann used the three- 
dimensional density profile suggested by the cosmological 
simulations of Navarro, Frenk & White (1996), namely: 



p(r) = 



r(r c + r) 2 



(3.21) 



where r is the three-dimensional radius and p s and r c are 
density and length scales. He projected this to obtain the 
surface density and hence the lens model. He showed that if 
both the positions and redshifts are available for the tangen- 
tial and radial arcs, then it is possible to find values for the 
density and length scales. If the redshifts of the arcs are not 



known, it is possible to constrain them. In particular, Bartel- 
mann (1996) pointed out a problem with the density profile 
(3.21) - for tangential arcs lying within the scale radius, the 
lens model predicts large values for the radial magnification, 
which in turn implies that the corresponding sources must 
be surprisingly thin in the radial direction. 

Can the more general family of surface densities (3.1) 
alleviate the problem described by Bartelmann? In addition 
to the arbitrary density and length scales (in our notation 
M and £ c ), we have the extra parameter 7 describing the 
cusp. By varying the value of 7, can we avoid the problem of 
large radial magnification at the tangential arcs? Equations 
(3.13) and (3.16) give the positions of the radial and tangen- 
tial critical curves for the family of lenses with convergence 
(3.11). Assuming that the sources producing the radial and 
tangential arcs are at different redshifts, we rewrite (3.13) 
and (3.16) as 



R 2 (l+R 2 -l=2M t R 2 -\ 
(1 + R 2 -~<) 2 =2M r ((l - 7R- 7 - R^), 



(3.22) 



where the subscripts on M indicate its redshift dependence. 
Defining the ratio a by 



Mt 

M r ' 



we obtain, using (3.22), 



<j{a) 



(l + R 2 t -")(l-l- 



Rf-'a 2 -^) 



o7(l 



-'Rt 



(3.23) 



(3.24) 



where a = R T /R t . The value of a is known from the posi- 
tions of the radial and tangential arcs. If the redshifts of the 
arcs are also available, then we can use (3.23) to find the 
corresponding value of a. Choosing a value for 7, we can 
solve (3.24) for R t . Since R t = £t/£c, and we can observe 
the value of the physical position of the tangential arc £ t , wc 
can obtain the scale £ c . Finally, (3.22) can be used to obtain 
the value of the density scale. 

Often the redshifts are unavailable and we must proceed 
more indirectly. As an example, let us consider the lensing 
cluster MS 2137, which is at a redshift of z c — 0.315. In this 
cluster, a tangential arc has been observed at £ t = 15". 5 
from the cluster centre. Fort et al. (1992) detected a radial 
arc in the same cluster at £ r = 5".0 from the cluster cen- 
tre. To proceed, we assume that £ c ~ 250/jT 1 kpc, a value 
obtained from numerical simulations of dark matter halos 
within the CDM cosmogony (cf. Bartelmann 1996). From 
the observed values of £ t and £ r , we conclude that a = 0.32 
and we find that Rt ~ 0.2. In Fig. 4(a), we have plotted a 
as a function of 7 for R t = 0.2 with a = 0.32 (and also for 
0.75a and 1.25a to give a feel for the likely uncertainties). 
By varying the value of 7, we can obtain a wide range in 
a. In Fig. 4 (b), we show a versus the redshift « r of the ra- 
dial arc for a number of different values of the redshift z t 
of the tangential arc. This plot shows that if a < 0.6, the 
redshifts of both sources must be unreasonably large. We 
should therefore choose a value of 7 for which a ~ 1. 

Another constraint on our choice of 7 is the radial mag- 
nification produced at the position of the tangential arc. 
Assuming the same value for £ c as before, we find that 
R T — 0.0645. Substituting this into (3.22) we can find M r as 
a function of 7 only. The value of the radial eigenvalue as a 
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Figure 4. (a) a versus 7 for different values of a. (b) a versus radial source redshift z r for various values of tangential source rcdshift 
z\. (c) Radial magnification at the location of the tangential arc plotted versus 7 for Rt = 0.2 and R r = 0.0645. 



function of position is then given by 

2M r ( 7 )((l- 7 ) J R r - 7 - J R?- 27 ) 



A r (i? r ) = 1 - 



(3.25) 



The radial magnification is just 1/A r . In Fig. 4 (c), the radial 
magnification at the position of the tangential arc is plotted 
as a function of 7. In order to obtain a value of the magnifi- 
cation which is of order unity, we must choose 7 in the range 
0-6 < 7 < 0.95. Fig. 4 (a) shows that values of 7 in the range 
0-6 < 7 < 0.7 correspond to acceptable values of a. In other 
words, the problem pointed out by Bartelmann is assuaged 
by a choice of a cluster lens that is more singular than the 
dark halo profile of Navarro, Frenk & White (1996). 



4 THE ISOTHERMAL DOUBLE POWER-LAW 
MODEL 

One of the double power-law models has a remarkable 
property - the lens equation can be re-cast as a solvable, 
quadratic equation for the image positions in terms of the 
source position. This is the model with an isothermal cusp. 
This property of invertibility is exceedingly scarce. It is pos- 
sessed by only two other known models - the Schwarszchild 
lens (Einstein 1936; Refsdal 1964) and the infinite isother- 
mal sphere (Kovner 1987a). 



4.1 The Image Positions 

When 7 = 1, the convergence diverges like R~ at small 
radii and like R~ 3 at large radii, viz; 



k(R) = 



M 



R(l + R) 2 



ip = 2M\og{l + R). 



(4.1) 



Deprojecting the surface density by the well-known Abel 
inversion formula (e.g., Binney & Tremaine 1987), we find 
that the spherical cluster has three-dimensional density 



p(r) = 



M 

2tt 2 
M 
2tt 2 



2r z + 1 



3arcosh(l/r) 



(1 



2r l 



(1 - r 2 ) 5 / 2 
3arcos(l/r) 



(r 2 



(r 2 - l) 5 /2 



r < 1, 



r > 1. 



(4.2) 





Figure 5. The variation of the time delay r and flux ratio of 
the images with separation of the source y from exact alignment 
with the isothermal double power-law lens. The curves are shown 
for total mass M = 1, 2, 3 and 4 respectively. As y — > 2M, the 
second image becomes fainter and fainter. It finally disappears 
when the source crosses the pseudo-caustic at y = 2M. At this 
point, the notion of a time delay loses its meaning. (Units used 
with £ 2 D S (1 + z d )/(cD d D ds ) = 1.) 



Here, r is the three-dimensional radius in units of the length 
scale. A careful Taylor expansion shows the density is regular 
at r = 1 and has the value: 



m M 



(4.3) 



The three-dimensional density diverges like r~ at small 
radii and like r~ 4 at large radii. The model is therefore akin 
to the famous halo model of Jaffe (1983). We shall refer to 
it as the isothermal double power-law model. The deflection 
angle at a displacement x is 



9 = w ■■ 



2Msign(x) 
1 + H ' 



(4.4) 



The circle \y\ = 2M is the pseudo-caustic. Source positions 
with \y\ < 2M give two images whereas those with \y\ > 2M 
give only one. The additional image is created or destroyed 
at the cusp. The pseudo-critical curves have degenerated 
into a point at the centre. 
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Table 2. Image Positions 
as a Function of Source 
Positions as Labelled on 
Figs 6 and 7. Roman nu- 
merals are used to indicate 
the quadrant in which the 
image lies. 



Si 


Lil 


s 2 


L 2 I, L 3 III 


S3 


L 2 I, L2lI,L3lI,L 3 III 


S 4 


LiIV 


s 5 


L 2 III, L 2 IV 


s 6 


LiIII, LiIV, L3III 


s 7 


L 2 IV, L 2 III 


s 8 


2L 2 III, L 2 IV, L4III 


s 9 


LiIII 


Sio 


L2III, L3III 



(e). 



7- L. 



1 2 

x, 



2^ S 4 



(f) 



1 2 



x 1 




I 


(b)_ 










i 





1 2 

7i 



Figure 6. The critical curve and caustic topologies for the 
isothermal double power-law model perturbed with external 
shear. In each case, the lens plane is shown on the left, the source 
plane on the right. In panels (a) and (b), C'i > C'2 > and there 
is one C point on each axis; (c) and (d), C\ > 0, > C'2 > —0.5 
and there is one C point on the X2-axis and one S point at the 
origin; in (e) and (f) C\ > 0, C2 < —0.5 and there is a C point 
and an F point on the X2-axis. The caustics and critical curves are 
drawn in full lines, the pseudo-caustics and pseudo-critical curves 
are drawn in broken lines. (Units used with M = 1.) 



X 1 



\ 


(c) 








L, 







1 

x, 




Taking y > without loss of generality, the images are 
located on either side of the source at 



x=\(y^l)±\A ± {y), 



with 



A ± (y) = y/(y±l)*+8M, 



(4.5) 



(4.6) 



provided the source lies within the pseudo-caustic. The im- 
age at x > is of type I, the image at x < is of type II. 
If the source lies outside the pseudo-caustic, then only the 
former image is present. When \y\ > 2M, the total magnifi- 
cation is 



"=2 



1 1 y(y+l)+4M 



2y A+(y) 



(4.7) 



Figure 7. The critical curve and caustic topologies for the 
isothermal double power-law model perturbed with external shear 
(continued). In panels (a) and (b), > Ci > —0.5 > C 2 and there 
is an F point on the 22-axis and an S point at the origin; in (c) 
and (d), —0.5 > C\ > C2 and there is one F point on each 
axis. The caustics and critical curves are drawn in full lines, the 
pseudo-caustics and pseudo-critical curves are drawn in broken 
lines. (Units used with M = 1.) 



whereas when Iwl < 2M, it is 



^Ty 



y(y + l)+4M y(y - 1) + 4M 



(4.8) 



A+(y) A-(y) 
The time delay r is also explicit. For any source position 
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and so scales roughly like M 1//2 . 



0.5 
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0.5 
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0.5 



(e) 




0.5 



Figure 8. As Fig. 6 but for an isothermal double power-law model 
with tiny core. This causes the replacement of the pseudo-caustics 
with true caustics and the destruction of the pseudo-critical 
curves. In each case, the central portion of the lens plane is shown 
on the left, the source plane on the right. In panels (a) and (b), 
C\ > C2 > and there is one C and one F point on each axis. 
Only the tiny inner radial critical curve is visible in the plot. The 
outer tangential critical curve is virtually unchanged from Fig. 6. 
In panels (c) and (d), Ci > 0,0 > C2 > —0.5 and there is one C 
point on the X2-axis and one F point on the xi-axis. In panels (e) 
and (f), Ci > 0, C2 < —0.5 and there is one F point on the x\ 
axis and a C point and two F points on the oi2-axis. (Units used 
with M = 1.) 



which generates two images, it is 
_gD s {l + z d ) 



D d D ds 
-y + 2M log 



±(W + l)A + (j,) + ±(y-l)A_(y) 
l+y + A+{y) 



l-y + A-(y) 



(4.9) 



where z d is the redshift of the lens. Fig. 5 shows the varia- 
tion of the time delay and the flux ratio of the images with 
separation of the source y from exact alignment. For the im- 
ages to be of comparable brightness, we require y < 1. The 
time delay is then approximately 



/2M$gj> B (l + Zd) 
D d D ds 



(4.10) 



4.2 Caustic Topology of the Perturbed Lens 

The aim of this section is to classify the structure of the 
caustics and critical lines of the perturbed isothermal model 
exactly. This job has already been performed for the Plum- 
mer model (the instance 7 = in (3.11)) in chapter 8 of 
Schneider, Ehlers & Falco (1992). The isothermal double 
power-law model has a pseudo-caustic not present in the 
Plummer model. The slightest breakage of the circular sym- 
metry causes the pseudo-critical curves to be displaced from 
the density singularity. These unusual features mean that it 
is well-worth classifying the caustic topology for this model. 

The lens action is considered as taking place in the pres- 
ence of an external larger-scale gravitational field with local 
surface mass density ^Ti + |F2 and local shear ^Fi — |r 2 . 
By choosing the orientation of the coordinates to diagonalise 
the tidal field, the lens equation becomes 



V = x 



1 - 



2M 



Fi 







r 2 



(4.11) 



Even in the presence of shear, the circle \y\ — 2M remains 
as the pseudo-caustic. The pseudo-critical curves are now 

2 _ 4M(i + r 2 ) - r(i + r 2 [r 2 - 2]) 



(r 1 -r 2 )(r[r 1 + r 2 -2]-4M)' 

where (r, 8) are polar coordinates in the lens plane. 
Let us introduce the shorthand 
M M 



Ci = 



■IV 



c 2 



(4.12) 



(4.13) 



Considering the first quadrant (i.e., x\ > 0, X2 > 0), critical 
points (C points) on the zi-axis occur at 



x x = i(l + 8C 2 )2 - \. (4.14) 

These exist provided C2 > 0. Folds (F points) occur at 

x 1 = yf = 2Cl-l, (4.15) 

provided Ci < — \. Alternatively, the fold can occur at the 
density singularity itself, when we refer to it as an S point. 
The discussion for the £ 2 -axis is similar with trivial changes 
of indices. 

Figures 6 and 7 show the topologies of both the lens 
(xi,X2) and source (2/1,2/2) planes for the distinct image con- 
figurations. In Figs 6 (a) and (b), there is one C point on 
each axis. In Table 2, the image regions corresponding to 
the different regions of the source plane are listed. For ex- 
ample, a source in region S3 produces four images, one in 
L2 in the first quadrant (L2I), one in L2 in the second quad- 
rant (L2II) and two in L3 in the second and third quadrants 
(L3II, L3III). Moving the source from S3 to S2 causes the two 
images in the second quadrant to fuse and only two images 
are left, one in the first quadrant (L2I) and one in the third 
(L3III). If the source is now moved across the pseudo-caustic 
from S2 to Si, the image in the third quadrant vanishes 
into the central density singularity, while the image in the 
first quadrant (Lil) remains. As expected, when the source 
crosses the pseudo-caustic, the number of images changes by 
one. Regions Li, L2 and L3 correspond to images of types 
1,1 and II respectively. There is no change in parity on cross- 
ing the pseudo-caustic. In Figs 6 (c) and (d), there is a C 
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point on the *2-axis and an S point at the origin. When- 
ever a critical curve ends at an S point, the corresponding 
caustic is discontinuous. As is evident in the diagram, the 
caustic terminates abruptly in the source plane as it touches 
the pseudo-caustic. Sources in S6 generate three images (see 
Table 2). One image is lost on crossing the pseudo-caustic 
from S6 to S5, whereas one image is gained on crossing the 
pseudo-caustic from S4 to S5. Regions Li, L2 and L3 cor- 
respond to images of types II, II and III respectively. In 
Figs 6 (e) and (f), there is both a C point and an F point 
on the 22-axis. Sources in Ss have four images (see Table 2). 
Traversing the pseudo-caustic from S6 to Ss or from S4 to S7 
causes the gain of a single image. Regions Li,L,2,L3 and L4 
correspond to images of types II,II,III and III respectively. 
If > Ci > C2 > —0.5 there are no critical points on the 
axes. Sources with \y\ > 2M possess one image, source with 
\y\ < 2M possess no images. The pseudo-critical curve lies 
at the origin. 

The remaining two possible topologies are shown in 
Fig. 7. They are included for completeness, but are not phys- 
ical, as there are regions of the source plane that correspond 
to no images in the lens plane. At first sight, this seems to 
contravene the theorem of Schneider (1984). If the time de- 
lay increases quadratically far from the lens, it seems that 
there must be at least one minimum of the Fermat surface. 
The reason why this theorem no longer holds here is that 
the values of the shear Ti and T2 now exceed unity. The 
Fermat potential decreases quadratically far from the lens 
and so there must be a highest point on the Fermat sur- 
face - but this corresponds to the density spike and not an 
image. A somewhat similar occurrence takes place in the 
singular point mass with shear or Chang-Refsdal lens (see 
Schneider, Ehlers & Falco, 1992, p. 261). In Figs 7 (a) and 
(b), there is an F point on the xi-axis and an S point at 
the origin. Figs 7 (c) and (d) show a similar configuration, 
but with an F point on each axis. In both cases, sources in 
S9 give rise to one image of type III in region Li. Cross- 
ing the pseudo-caustic from S9 to S10 generates a second 
image of type II in L3, while the pre-exisiting image moves 
from Li to L2. However, as the source crosses the caustic 
to Sn, both images merge and disappear into the central 
density singularity. Note that in Fig. 7 (b), the crossing of 
the pseudo-caustic from S9 to Sn is accompanied by the de- 
struction of a single image whereas the crossing from Sg to 
S10 is accompanied by the creation of a single image. Sources 
in Sn have no images at all. 

As an aid to understanding, it is helpful to introduce 
a tiny core into the isothermal double power-law model so 
that it becomes 



M 



(R+e)(l + R) 2 ' 
with corresponding deflection angle 
2Msign(x) 2Me 



e 



(l-e)(l + |x|) |*|(1 -e) s 



log 



6(1 + |Z|) 



(4.16) 



(4.17) 



The changes in the topology of the caustics and the critical 
curves are shown in detail in Fig 8. This figure is directly 
comparable with the cases examined in Fig 6. The introduc- 
tion of the infinitesimal core causes the conversion of the 
pseudo-caustics into true caustics, as is evident on compar- 
ing Figs 8 (b) and (f) with Figs 6 (b) and (f). Now, the 



number of images changes by two on crossing these curves. 
The pseudo-critical curves are completely destroyed by the 
introduction of the core. Replacing them are tiny critical 
curves generated near the origin. For example, Figs 8 (a) 
and (e) show the emergence of a minute radial critical curve 
from the origin because of the infinitesimal core. Only the 
inner portion of the lens plane is shown in the figures - the 
outer critical curve is virtually unchanged from the singular 
model in Figs 6 and 7. So, not all the critical points on axis 
are visible in the details of the central regions. In Fig. 8 (d), 
only part of the pseudo-caustic is converted to a caustic, 
the remaining part being destroyed. In the lens plane, the S 
point of Fig. 6 (c) becomes an F point in Fig. 8 (c). 



5 THE POWER-LAW LENSES 

In this section, we study density cusps of infinite extent with 
elliptic equipotentials. Lenses with equipotentials stratified 
on similar concentric ellipses have been studied before (e.g., 
Kovncr 1987b; Blandford & Kochanek 1987; Kochanek & 
Blandford 1987; Kassiola & Kovner 1993; Witt 1996; Witt 

6 Mao 1997), although the emphasis has generally been on 
models with softened cores. By contrast, our interest is fo- 
cused on the lenses with central density singularities. 

5.1 The Potential-Convergence Pair 



The deflection potential is 



(5.1) 



where q is the axis ratio of the equipotentials, 7 is the cusp 
index and the constant A fixes the overall scaling. All the 
models have infinite total projected mass. 

As Contopoulos (1954) noted, three-dimensional distri- 
butions that are stratified on similar concentric ellipsoids 
project to two-dimensional distributions stratified on similar 
concentric ellipses. The projected potential (5.1) therefore 
corresponds to the three-dimensional galaxy models written 
down by Evans (1993, 1994) and called the power-law galax- 
ies. Schneider, Ehlers & Falco (1992) caution that models 
with elliptical equipotentials may not be very realistic, so it 
is worth pointing out that the power-law galaxies have been 
successfully used to represent the nearby elliptical M32 (van 
dcr Marel et al. 1994) as well as the inner 500 parsecs of the 
Galactic bulge (Evans & de Zeeuw 1994). The convergence 
is cusped like 



A(2- 7 ) 



ax\ + bx\ 



2q 2 (x? + ^-2)1+7/2 



with 



7), b 



2 (l-7)- 



(5.2) 



(5.3) 



The axis ratio of the equipotentials q describes the flattening 
of the model and can be restricted to lie in the range < 
q < 1 without loss of generality. The projected mass density 
of the singular models has ellipticity e 

„2 



7+1 1 1/7 



l + ? 2 (l-7) 



(5.4) 



The convergence is positive if 7 < 1 + q 2 (e.g., Blandford 
& Kochanek 1987). If the model is viewed as a projected 
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Figure 9. The critical values of q for naked cusps (q n ) and the 
umbilic catastrophe (q u ) plotted as functions of 7 

power-law galaxy, the constraints that the three-dimensional 
density and that the two-integral distribution function are 
positive definite are more severe and given in Evans (1994). 
The parameter 7 describes the rotation curve of the model. 
If 7 = 1, the model has an asymptotically fiat rotation curve 
(e.g., Binney & Tremaine 1987; Evans 1993). If 7 > 1, the 
rotation curve is declining, whereas if 7 < 1, the rotation 
curve is rising. The range of 7 is normally restricted to < 
7 < 2. The lower limit comes from requiring the convergence 
to vanish at spatial infinity, the upper limit is the Keplerian 
case of a point mass. 

The tangential critical curve has intercepts 



2 

*i,t 



Aq- 2 {2- 1 ) 



2/ 7 



q 2 A(2 - 7) 



2/7 



(5.5) 



(5.6) 



The condition for the existence of a tangential critical curve 
is 



,4(2 -7) > 0. 



(5.7) 



It is always satisfied. If the cusp is strong (7 > 1), there is a 
transition from two to four images on crossing the tangential 
caustic. If the cusp is weak (7 < 1), there is a transition from 
one to three images. The radial critical curve has intercepts 



2 



*2,i 



A(2-7)(l- 7 ) 



2/7 



q 2 AT 2 (2- 7 )(l- 7 ) 



2/7 



(5.8) 



(5.9) 



The radial critical curve exists only if 7 < 1. Quintuple 
imaging occurs if 



A(2-7)(l- 7 )>0. 



(5.10) 



That is, quintuple imaging can occur for all power-law galax- 
ies with weak cusps. The critical curves become smaller both 
in the limit of increasing homogeneity (7 — > 0) and in the 
Keplerian point mass limit (7 — > 2). 

In the source plane, the semi-axes of the tangential and 
radial caustics are 



3/2,t 

3/1 ,t 



3/l,r 
3/2,r 



= Q 



2/7-1 



(5.11) 



For 7 < 1, there are two important topology changes that 
can occur. For models only weakly perturbed from circular 



symmetry, the radial caustic lies entirely outside the tan- 
gential caustic. Models flatter than the critical flattening q n 
given by 



2 1 

<7n = 1 



4n /7 7(l - 7) 



1/7-1 



(5.12) 



possess "naked cusps" . That is, two of the cusps of the tan- 
gential caustic lie outside the radial caustic. The umbilic 
catastrophe occurs at a critical flattening q u given by 



2 
<7u 



■ 7- 



(5.13) 



Fig. 9 shows the behaviour of the critical axis ratio q n and 
g u for the onset of naked cusps and the umbilic catastrophe 
respectively as a function of cusp index 7. 

For the moment, let us consider the model without ex- 
ternal shear. In general, the lens equation must be solved nu- 
merically. However, when the source is on-axis and aligned 
with the lens, the image positions and magnifications can 
be found exactly. This is sufficiently uncommon for non- 
axisymmetric lenses that it is worth giving the solution ex- 
plicitly. There are always two images on the minor axis at 
(0, ±£2), where 



x 2 — q 



A(2- 7 ) 



1/7 



These have even parity and the magnification fx is 
1 

A* : 



7(1 - q *y 



(5.14) 



(5.15) 



There are always two images on the major axis at (±xi,Q), 
where 



Xl = A{2 - 7) 



1/7 



These have odd parity and the magnification fi is 
a 2 

A* = 



7(1 -q 2 )' 



(5.16) 



(5.17) 



Only these four images occur when the cusp is strong. For 
weak cusps, there is an additional (infinitely de-magnified) 
image of even parity at the origin. The total magnification 
Atot is 



2(1 + g 2 ) 
7(1 -q 2 Y 



(5.18) 



This equation is exact when the source is on-axis and is ap- 
proximate when the source lies within the tangential caus- 
tic. This equation has previously been given by Kassiola & 
Kovner (1993) for the special case 7 = 1. 

The 7 = 1 model has the characteristic properties of 
an isothermal cusp. Unlike the perturbed circular model of 
Section 4.2, it has a truly non-axisymmetric pseudo-caustic. 
Its equation can be calculated by the methods of Section 2.2. 
If the phase of z = x\+ix 2 is 6 and the phase of z s = y\ -M3/2 
is 4>, then the equation of the pseudo-caustic is 



1 1 2 2.2 
W =3/1+3/2 



\OXiJ \OX2J 

A 2 (cos 2 e + q- 4 sm 2 6) 



(5.19) 



> + q- 



As the source crosses the pseudo-caustic, an extra image 
of type II is generated at the density singularity. On the 
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Figure 10. The critical curve and caustic topologies for the 
power-law lenses in the presence of external shear. In each case, 
the lens plane is shown on the left, the source plane on the right. 
This figure is for the case of a strong cusp with 7 = 1.25 and 
q = 0.6. [Figs 10 (a) and (b) have Ti = 0.1, T 2 = 0.5, Figs 10 (c) 
and (d) have IT = 0.8, T 2 = 0.7] 

pseudo-caustic, we know from the lens equation that the 
relation between the two phases is 



tan 6 = q tan <j>, 

and so the pseudo-caustic is an ellipse of form 



2 . 22 



(5.20) 



(5.21) 



The pseudo-caustic is concentric with the elliptic equipoten- 
tials, but it is elongated in the opposite direction. For this 
model, the lens equation can be reduced to a quartic equa- 
tion. The quartic degenerates to simpler quadratics if the 
source lies either on the major axis or on the minor axis. 
Suppose the source lies on the major axis at (3/1, 0), with 
2/i > without loss of generality. Then there is one image if 
y i > A at 

X!=yi + A, x 2 = 0. (5.22) 
Initially, this image has even parity. Its magnification is 
A + yi 



/' : 



\ yi -A(q-*-l) 



(5.23) 



On crossing the pseudo-caustic at yi = A, there is a second 
image generated at 



xi = yi — A, X2 — 0. 
It has odd parity and finite magnification 
A-yi 



A(q~ 2 — 1) + 1/1 



(5.24) 



(5.25) 



On crossing the tangential caustic at yi = A(q 2 — 1), the 
image (5.22) changes from even to odd parity. Simultane- 
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Figure 11. As in Figure 10, but for the case of a weakly 
cusped power-law lens with 7 = 0.75 and q = 0.6. [Figs 11 
(a) and (b) have IT = -0.5, T 2 = -0.1, Figs 11 (c) and 
(d) have F 1 = 0.3, T 2 = 0.1 and Figs 11 (c) and (f) have 
IT = -0.7, r 2 = 0.1] 



ously, two new images of even parity are generated at 



Xl 



yi 



X 2 



(l-q 2 Y q 
The magnification of each image is 

A 2 q- 2 (q- 2 - 1) 



{q- 2 - I) 2 



1/2 



A 2 (q~ 2 - l) 2 



(5.26) 



(5.27) 



Note that a change in parity of one of the pre-existing im- 
ages must occur on crossing a pseudo-caustic for our general 
theorems of Section 2.2 to hold good. The images and mag- 
nifications when the source lies on the minor axis can be 
found similarly. 

5.2 Caustic Topology of the Perturbed Lens 

Figures 10 and 11 show the possible topologies of the caus- 
tics for the power-law lenses in the presence of external 
shear. In this case, the lens equations become 

A(2- 7 ) 



yi = xi 



l+Ti- 



x 2 + x\q~ 



2)7/2 



+ T2X2, 



(5.28) 
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Table 3. Image Po- 
sitions as Function 
of Source Positions 
as Labelled on Figs 
10 and 11. 



Si 


Li , L2 


s 2 


2Li, 2L 2 


s 3 


2L 3 


s 4 


3L3, L4 


s 5 


L 5 


Se 


L5, L6, L7 


s 7 


2L5, 2Lq, L7 


Ss 


2L5, L8 



and 

2/2 = x 2 



Aq- 2 (2- 1 ) 



+ r 2 xi, 



where the shear matrix is taken as traceless and of form 



(5.29) 



ri r 2 
r 2 -ri 



(5.30) 



These equations have been recently considered by Witt & 
Mao (1997), who prove a number of interesting theorems re- 
garding the image positions. Since the shear axes do not, in 
general, coincide with the major and minor axes of the lens 
potential, the symmetry axes of the caustics are also differ- 
ent from those of the lens. Consequently, Table 3 - which 
gives the image regions corresponding to different regions in 
the source plane - does not contain information about which 
quadrants will contain the images. The quadrant in which a 
given image is found varies according to the exact position 
of the source. 

Fig. 10 shows the two distinct topologies which occur 
for a strong cusp (7 = 1.25). When the external shear is 
small compared to the lens shear, there is only one critical 
curve in the lens plane and a corresponding caustic in the 
source plane, as shown in Figs 10 (a) and (b). A source in 
region S 2 produces four images, two in region Li of the lens 
plane and two in region L 2 . Moving the source from S 2 to 
Si causes two of these images to fuse together leaving only 
two images, one in region Li and the other in region L 2 . If 
Ti and T 2 satisfy 



r? +rl 



(5.31) 



then the caustic becomes extremely elongated along the axis 

r 2 



in 



(5.32) 



i + rr 

The result of further increasing the shear so that Y\ +T% > 1 
is shown in Figs 10 (c) and (d) - there are now two critical 
curves and two caustics. A source in region S4 produces four 
images, three in region L3 and one in region L4. If the source 
if moved from S4 to S3, two of the images fuse and only two 
images are left, both lying in region L3. It is also possible to 
have a situation in which there is no caustic. 

Fig. 11 shows the three possible critical curve and caus- 
tic topologies for the case of a weak cusp (7 = 0.75). Figs 11 



(a) and (b) show the radial caustic lying entirely outside the 
tangential caustic. Alternatively, as in Figs 11 (c) and (d), 
the tangential critical curve can completely contain the ra- 
dial caustic. For intermediate values of the shear, the tangen- 
tial curve can lie partly outside the radial curve and partly 
inside, as shown in Figs 11 (e) and (f). If we continue to in- 
crease the magnitude of the external shear, the outer curve 
in Fig. 11 (d) is observed to be stretched along its major 
axis while being simultaneously compressed along its minor 
axis. When Fi and F 2 satisfy (5.31), this curve becomes a 
straight line. Outside the unit circle in the (Fi,r 2 ) plane, 
there exists only one caustic, with two cusps. 

In Figs. 11 (a) and (b), a source placed in region S7 pro- 
duces five images as expected for a weak cusp. There are two 
in region L5, two in L6 and one in L7. If the source is moved 
across the caustic to region S6 , one pair of images fuses to- 
gether to leave three images, one in each of the regions L5, 
L6 and L7. Finally, if the source is moved to region S5 an- 
other pair of images fuses to leave a single image in region 
Lb. The locations of the images corresponding to the source 
regions in Figs. 11 (d) and (f) can be found from Table 3. 
As described above, if Ti and T 2 lie outside the unit circle 
in the (Fi,r 2 ) plane, the tangential critical curve no longer 
exists. Sources placed inside the one remaining caustic give 
three images (two lie outside the critical curve and one lies 
inside) while those placed outside give only one image lying 
outside the critical curve. 



6 CONCLUSIONS 

This paper has studied the properties of lenses with central 
density cusps. These are models in which the convergence k 
diverges as a function of projected radius R like 



k oc R 



(6.1) 



near the centre. Numerical simulations (e.g., Navarro, Frenk 
& White 1996), dynamical arguments (e.g., Evans & Collett 
1997) and observations (e.g., Lauer et al. 1995; Faber et al. 
1997) all suggest that galaxies and clusters commonly have 
such singular density profiles. By generalising the index the- 
orem to allow for central density singularities, we suggest the 
following classification of cusps on the basis of their lensing 
properties: 

(1) Strong density cusps (1 < 7 < 2) give rise to an even 
number of images. There are as many images of odd parity 
as there are images of even parity (m — mi + nm = 0). 
Models with strong cusps do not possess radial caustics. 

(2) Weak density cusps (0 < 7 < 1) give rise to an odd 
number of images. The images of even parity outnumber the 
images of odd parity by one (ni-nn+nm = 1). Models with 
weak cusps possess radial caustics. Even though such models 
are singular, they can produce radial arcs (cf. Mellier, Fort 
& Kncib 1993; Bartelmann 1996). 

(3) Isothermal density cusps (7 = 1) always possess pseudo- 
caustics. If the source is within the pseudo-caustic, there 
is an even number of images (m — «ii + n-n = 0). If the 
source is outside the pseudo-caustic, there is an odd number 
(ni — Tin + nn = 1). A method for computing the loca- 
tion of the pseudo-caustic is given. Note that the pseudo- 
caustic differs from a true caustic in two ways. First, the 
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magnification of a point source is finite, not infinite, at the 
pseudo-caustic. Second, the number of images created or 
destroyed on crossing the pseudo-caustic is one, not two. 
A circular lens, even if perturbed by external shear, has a 
pseudo-caustic that is a circle. This is not the case if the lens 
is truly non-axisymmetric. 

As specific examples, we have presented two families of 
models. The double power-law family has a central cusp, 
a transition region and an outer envelope. The power-law 
family are a set of elliptical cusps of infinite extent. They 
are projected power-law galaxies. In both cases, detailed 
properties of the individual lenses are given, together with 
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a classification of the caustic structures in the presence of 
external shear. As an application, we consider a problem 
pointed out by Bartelmann (1996). He suggested that lens- 
ing may be problematic for the cusped Navarro, Frenk & 
White (1996) density profile in cases where both radial and 
tangential arcs are present. On matching the positions of the 
arcs, this model predicts large values for the radial magnifi- 
cation at the tangential arc, which in turn implies that the 
corresponding sources must be radially very thin. We have 
shown that this problem is not so serious for lenses that 
have more singular density cusps - and, indeed, higher res- 
olution simulations of halo formation suggest that the cusp 
slope may have been originally underestimated (Fukushige 
& Makino 1997; Moore et al 1997). 

Our investigation has uncovered a new model - the 
isothermal double power-law model - for which the lens equa- 
tion is exactly solvable for any source position. This scarce 
property is only possessed by two other known circularly 
symmetric lenses, the Schwarzschild lens and the isothermal 
sphere. The isothermal double power-law model has a three- 
dimensional density which behaves like r~ 2 in the centre and 
like r -4 in the outer parts, so it is a realistic model for a 
galaxy or cluster of finite extent with a flat rotation curve. 
As the image positions are explicitly available, the lensing 
cross-sections are straightforward to calculate. 
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